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ABSTRACT 

— ». , We discuss the relevance of the curvature radiation of pions in strongly magnetized 

^^ I pulsars or magnetars, and their implications for the production of TeV energy neutrinos 

^^ ' detectable by cubic kilometer scale detectors, as well as high energy photons. 
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^S^ , 1. Introduction 

O ' In pulsar and magnetar polar caps one expects the formation of vacuum or space-charge limited 

-*^ ■ gaps (e.g. Harding and Lai, 2006), in which an electric field component parallel to the magnetic 

C^ . field accelerates particles to relativistic velocities. The particles move away from the caps along the 
magnetic field lines, producing photons by curvature radiation, which at some height h produce an 

K^ , e^ pair front that determines the upper extent of the gap (e.g. Ruderman and Sutherland, 1975, 

'!^ I Harding and Muslimov, 2002). The caps with 0,-B < can accelerate positively charged protons or 

—.^ ' ions, whose interaction with X-rays from the polar cap will lead to very high energy (VHE, £'> 10^^ 

Tij- ■ eV) neutrinos through photomeson interactions, p7 -^ vr"'" — > ^"'"f^ -^ e'^VeV^ (Zhang et al, 2003). 

[*-. ■ Here we explore the consequences for neutrino production of a different mechanism, the curvature 

o 

OO 

O 



radiation of pions by protons. 



Relativistic protons interacting with a strong magnetic field can produce pions, a quantum 
treatment of this process being given by Zharkov (1964). Ginzburg and Zharkov (1965) used a 
/\ ' semi-classical approach to emphasize the analogy between this pion radiation process and the usual 

C^ ■ synchrotron radiation of photons by protons. This semi-classical method was further explored 

for vr*^ and p mesons by Tokuhisa and Kajino (1999) in the context of pulsars. The pion radiation 
mechanism is characterized by a parameter x = i'^{B/Bq) similar to that used for electron processes 
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but with the proton mass replacing the electron mass, with Bq = m'i(? /eh = 1.5 x 10^*^ G. They 
showed that the synchrotron-like vr*^ radiation process becomes competitive with (proton) photon 
synchrotron radiation when the above parameter x^O.l. The related synchrotron-like vr"*" radiation 
by protons in a magnetic field requires, however, a different semi-classical treatment, since the 
protons become neutrons in the process, and this was calculated by Herpay and Patkos (2008). 

In the astrophysical literature, particle acceleration in pulsars or magnetars is considered both 
near the polar caps ("inner" or polar gap models), and further out near or beyond the light cylinder 
("outer" gap or wind models); e.g. see review by Harding and Lai (2006). Examples of the former 
are e.g. Blasi, Epstein and Olinto (2000), Arons (2003), etc.; while examples of the latter are, e.g. 
Sutherland (1979), Harding and Muslimov (2001, 2002, 2005), Thompson (2008), etc.. Here, we 
consider specifically the case of protons or ions accelerated in the inner regions, just above the polar 
caps of magnetars. In the inner acceleration region discussed here, the accelerated particles are 
expected to radiate away their transverse energy so efficiently that they are effectively constrained 
to the ground Landau level as they move along the curved field lines. The principal photon radiation 
loss is then curvature radiation (Ruderman and Sutherland, 1975; Harding and Lai, 2006). In this 
situation, the synchrotron-like pion radiation (which depends on the proton transverse energy) will 
be similarly suppressed. However, as for photon curvature radiation, the motion along the curved 
field lines produces an acceleration equivalent to a Larmor motion in a fictitious magnetic field, 
which will lead to curvature pion radiation (Berezinsky, Dolgov and Kachelriess, 1996). Unlike 
previous authors, in this paper we concentrate on vr^ production by the curvature-like mechanism 
of both protons and heavy nuclei, occurring in the immediate vicinity of the polar caps, and we 
explore its consequences for the production of VHE neutrinos, as well as their detectability by 
cubic kilometer Cherenkov detectors such as IceCube and KMSNeT. We compare the resulting 
neutrino fluxes to those from the photomeson process, and discuss the associated cosmic ray as well 
as GeV-TeV photon production. 

In ^we discuss an approximate acceleration and radiation loss model, including both photon 
and pion radiation losses. Simple time-scale arguments are presented to estimate the maximum 
particle energies possible and the ratio of photon to pion losses. In ^we consider in more detail the 
equations of motion in self-consistent polar cap magnetic field models, and calculate numerically 
the acceleration and energy loss rates in a general manner. These results are in good agreement 
with the approximate estimates of the previous section. In §3 and §4 we estimate the resulting 
neutrino luminosity from curvature pion decays and from photo-meson effects. In ^ we discuss 
the possibility of detecting these effects in kilometer scale neutrino detectors and at GeV photon 
energies. 



2. Ion Acceleration and Energy Losses 

The maximum potential drop in a magnetar of radius a = W^qq cm, rotation frequency 
Q = 2-k/P (where P = 10~^P_3 s is the period) and surface magnetic field at the pole B = 10^^i?i5 



G is <& = {Vt^Bo? /2c^) ~ (6.6 x 10^^ \)Bi^a\Pzl- However, the self-consistent gap structure of 
magnetars is not well known, e.g. for the Q ■ B < Q caps where proton acceleration can occur. 
Vacuum gaps may be expected (Medin and Lai, 2008) under some assumptions about the surface 
cohesive energy and thermal structure. On the other hand, the uncertainties and approximations 
involved leave open the possibility of space-charge limited flow (SCLF) gaps, and in this paper 
we assume that the latter situation applies. Ions may be lifted from the magnetar surface and 
accelerated by an unsaturated electric field component parallel to the magnetic field lines for SCLF 
gaps (e.g. Harding and Muslimov 2002) up to a height similar to the polar cap radius ~ r^c = 
(27ra^/cP)^'^ ~ 4.6 x lO^Og -P_3 cm. Beyond a height h ~ rpc, ions may further be accelerated 
by an electric field in the saturated regime, given by ^^y ~ 5 x lO^^BisPTg in c.g.s. units (see, e.g., 
eq. 1 in Harding and Muslimov 2002). For heights much smaller than the stellar radius, /i <C a for 
which the saturated field is constant, the energy e reached by the particle is e = ZeEuh, during an 
acceleration time 

ta = e/cZeE\\ ~ 2.2 x IQ-^egSf/P^^gZ"^ s, (1) 

where eg = (e/lO^eV). Charged particles accelerated from the polar cap move following closely 
the field lines, the main photon radiation loss being curvature radiation, whose power is Vc = 
(2e^c7^/3i?^). Curvature radiation generally dominates over (transverse) synchrotron radiation at 
low heights in the polar caps. For the limiting open field line of a neutron star rotating with period 
P = 10^^P_3 s defining a light cylinder Ri = cP/2Ti the curvature radius is 

R^ = (4/3) (ai?L) 1/2 ^ 3 ^ 10^P%\1^^ cm. (2) 

and the photon curvature radiation loss time for protons of energy e is 

tc = {2,eRl/2e^c-i'^) = 2.4 x lO^^eg^^P-sOe s. (3) 

Another energy loss mechanism for protons in a strong magnetic field is the synchrotron-like pion 
radiation losses (Ginzburg and Zharkov, 1964). In terms of the parameter x defined as 

X = i{B/Bq), (4) 

where Bq = mic' /eh = 1.5 x 10^*^ G, a semi-classical calculation in the x ^ 1 asymptotic limit 
gives the energy loss rate 'Ppy^ = dEp/dt of a proton due to tt^ emission as 

VAx^^)^{9V(!>)myh~\ (5) 

which in quantum mechanical calculation acquires a further x^ factor (Ginzburg and Zharkov, 
1964). Here g^ /he ~ 14 is the strong coupling constant. The corresponding low x limit is 

^7ro(x < 1) = (5'2/\/3)m^mpC^^~2xexp[-(\/3/x)(m^/mp)] (6) 

(Ginzburg and Zharkov, 1964; Tokuhisa and Kajino, 1999). For x,^0-l these losses can exceed 
those from photon synchrotron radiation. In the related process of tt~^ meson radiation, the proton 



becomes a neutron during the emission process, and a different semi-classical calculation method 
is needed, as done by Herpay and Patkos (2008). A quantum calculation of this process is that of 
Zharkov (1964). In the limit of x 'C 1 which is of interest for us here, the proton energy loss rate 
due to 7r+ radiation is 
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However, for the same reason that photon curvature radiation dominates over photon synchrotron 
radiation, the pion radiation due to the proton transverse momentum will be dominated by pion 
radiation due to the longitudinal component along the curved field lines (Berezinsky et al, 1996). 
We define an equivalent 'curvature" B±^ field perpendicular to the particle trajectory, which would 
give a radius of curvature Re equal to the gyro-radius of a proton of energy e and Lorentz factor 7, 



S_L = e/eRc = 1.1 X lO^egP: 



-1/2-1/2 p 



The corresponding 'curvature" x parameter is 

X± = i{B^/Bq) = ^\h/R,mpc) ~ 7.84 x IQ-^e^p 



-1/2 -1/2 



and in terms of this parameter the photon curvature power is 



V. 



7,c 



(2e2c7V3i?e) = {2e'^/3hc){my/h)x^ . 
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From now on we will refer to x = X±) s-nd in our case we will be interested mainly in the x *S 1 
limit. Within the range of this inequality it is reasonable to assume that the small momentum 
transfer regime applies and provides an adequate approximation to the interaction rate. This is the 
region where the simple effective models leading to equations ([6]) and ([7]) are valid. This intuitive 
description is based on the assumption that the transverse momentum of the emitted pions is of 
the same order of magnitude as that of the proton moving along the curved trajectory in the strong 
magnetic field (Berezinsky et al, 1996). An accurate discussion of this issue is still an open question. 
In this X *^ 1 limit, the cooling time against pion curvature radiation of a particle of energy e, 
tTz,c = e/T^TT.c, using eq. ^ with eq.®, is 
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[4.27 X 10-^^€g'P_'^aQ' exp[4M6€g^Py^a, 
where eg = (e/lO^GeV) is the particle (longitudinal) energy. 
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(11) 



For the case of heavy ions (A,Z), e.g. Fe (56,26), for an ion of energy e the ion Lorenzt 
factor is 7^ = e/{AmpC^), and each constituent proton moves with 7^. We can use the previous 
equations for the photon curvature losses 7^7,0 and in principle also the curvature pion losses Vt^^c 
of the individual constituent protons moving with 7^, and find the corresponding losses for the 



ion as 7^7, c,A = ZV^^ci'fA)- The proportionality factor Z assumes incoherent radiation by the 
constituent nucleons, since at high energies the dominant wavelength becomes smaller than the 
average separation of the constituents. Thus 'P-y,c,A = ZV^^c = Z{2/2>){e'^ /hc){mic^/K)x^. For an 
ion energy e^ = 10^^ en GeV, the acceleration and incoherent photon curvature times are 

ta,A = 8.6 X lO-^ii^fsiP-s^aV s> (12) 

i^,,,A = 1.07 X 10-9eri'^_3a646^26' S' (13) 

where we used the same curvature radius Re of equation ([2|) as for protons, since the nuclei move 
along the same field lines. The fictitious 'curvature" magnetic field B± and x± for ions is the same 
eq. ([8]) and eq. ([9]) as for protons, except for using ^a-, 

XA = 1a{BJBq) = 2.62 X lQ~\l^Al^{a(,P_:,)-^l^ . (14) 

If the ions did not fragment before reaching the appropriate energy threshold, the pion curvature 
radiation of heavy ions {A, Z) moving along the field lines would, in principle, be due to the 
individual nucleons in the nucleus moving with Lorentz factor ^a and with the x parameter of eq. 
(fnil . The pion emission power would then be Vt^^^a = ^"Ptt.c for the above XA-, the protons in 
the nucleus radiating tt"*" and the neutrons radiating vr^. However, as discussed in §3.1 and §3.2 
, the interaction with the curved magnetic field and the radiation recoil of individual nucleons is 
found to exceed the binding energy of the nucleus, leading to its fragmentation, before reaching the 
energy at which the nucleus could radiate pions. The individual resulting protons would, however, 
continue to be accelerated and would radiate individually, as described before. 

For protons, it is seen from Fig. [T]that the intersection of the vr" energy loss timescale with the 
acceleration time occurs generally at higher energies than the intersection of the tt"*" losses and the 
acceleration. The shorter tt"^ time is mainly due to the lower energy cutoff in the radiation power 
Vt^+ caused by an {111^+ /rrip)^ factor in the exponent, versus an m^o/mp in the exponent of the vr*' 
power. Thus, the proton energy is determined by a balance between acceleration energy gains and 
a combination of the 7r+ radiation and photon curvature losses. Assuming that the particles can 
reach the corresponding height h (which is plausible in the saturated field regime r^^ ■^h -^ a), the 
maximum particle energy for a given value of B and P is approximately given by the intersection 
of the appropriate acceleration time with either the 7r+ or the curvature loss time curve, whichever 
occurs at the lowest energy. Thus, from Fig. [H for B = 10^^ G and P = 10~^ s the curvature 
pion losses are beginning to get competitive with curvature photon losses, with t^+/tc ~ 0.2 at 
(-max ~ 0.65 X 10^ GeV. The pion losses dominate for Bi^PZ-^ >2. However for lower fields or longer 
periods -BisPjo <,1, it is the photon curvature losses that determine emax- 



3. Trajectories, Acceleration and Energy Losses from the Solution of the 

Proton/Ion Equations of Motion 

In the previous section, we determined the maximum energy of the protons/ions by comparing 
the characteristic timescales of different radiations with that of the acceleration calculated without 
taking into account the back-reaction of the radiation. Moreover, in estimating the characteristic 
times we made the further assumption that the perpendicular momentum of the accelerated particle 
to the magnetic field lines is much smaller than the parallel component. The simplicity of this 
physical picture greatly helps our intuition, still its fundation needs some elaboration. Therefore it 
is interesting to compute the trajectory and the saturation energy of the charged particles directly 
from the full equation of motion which also takes into account the backreaction of the radiation. 
In this section we discuss in detail the proton/ion trajectory through the acceleration region near 
the polar cap. 



3.1. Equations 

In a given electromagnetic field, the relativistic equation of motion of a proton is 



7 = — ^+^T=f7X^)-^^^7-^^ 7=^> (15) 
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where we added on the right hand side the last term to the usual Lorentz force to take into account 
the energy loss of the particle {V = dE/dt) due to radiation processes (photon, vr*^, vr"^). The 
direction of this back reaction force is chosen opposite to its momentum, which is true for large 
I7I (a similar description of the back reaction of the radiation was proposed in Takata et al (2008) 
and Harding, Usov and Muslimov (2005)). Here, V results from the curvature of the trajectory. 
It is not split up into separate synchrotron and curvature (along B) radiation parts corresponding 
to transverse and parallel momentum components, respectively. This approach enables us a more 
complete description, because the presence of an accelerating electric field continously modifies the 
curvature even if E\\B. Therefore in eq. (fT5]) V depends on the total curvature radius (Re), which 
defines a curvature parameter x (Berezinsky et al, 1995), 

X = -^-—i (16) 
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where Re can be expressed with 7 and 7, 



1 _ 1 vi + 7^ ^ ._ivr+f2/— 



7x7 



c I7P 



1'r - {iif- (17) 



Re c I7P 
Here P is the sum of the power of 7 and tt radiations given in section 2: 

V = Vc + V^o+V^+=v{\%x{l,1)), (18) 



The formulae are applicable only for x < m.,^/mp. It can be seen from the above expressions that 7 
appears on the right hand side of the equation of motion ()15p in a very complicated way. However if 
the ratio of the power of the energy loss to the energy is small, one can solve the equation iteratively, 
e.g. one substitutes 7 into x with its value calculated from the equation of motion with only the 
Lorentz term (without V). 

In this approximation one considers six coupled first order differential equations for r and 7 
wherein time derivatives appear only on the left hand sides 



"^ 7 = — (m + ^^7 X m] - ^(1^1 '^(J'^")) v^ ^, (19) 
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One can easily check that this approximation gives back the usual formulas for the intensity of 
the synchrotron or curvature radiation as long as E = or B = 0, respectively. Note that 
equations in (J19p would be further simplified if £^ 1 1 i? (see Takata et al, 2008) using appropriate 
local coordinates to treat the motion parallel and perpendicular to the magnetic field lines. In the 
general case {E x B ^ 0) one cannot avoid to calculate all three independent components of 7. 
Then the use of a coordinate system with one axis fixed along the magnetic field lines does not 
bring any extra advantage over the use of a unique coordinate system for the whole motion. 

Unfortunately the exact expression of the electromagnetic field to be used in (|15p is not known 
fully since the field-charge, the motion and emission of any charged particle modifies the accelerating 
field. The full selfconsistent description of the inner gap is very complicated as was already shortly 
discussed in section 2. In our actual calculations we considered a dipole field for the magnetic 
part. In addition a so-called unsaturated expression was used for the electric field close to the 
surface of the magnetar and saturated field form is adjusted to it far away from the surface. The 
explicit expressions of the electric field strength in the two regimes are given by Harding and 
Muslimov (2001, 2002). These electric fields take into account the infiuence of the accelerated 
electrons/positrons and their pair production, which may be the most important effect for the 
formation of a self consistent electric field. The corresponding unsaturated and saturated potentials 
can be found in equations (20) and (24) of Harding and Muslimov (2001). 



In the course of the numerical solutions of ()19p we switched from the unsaturated to the 
saturated field at the altitude where the component of the unsaturated electric field parallel to the 
magnetic field equals to the corresponding component of the saturated field. 

By consulting the literature it appears, that our choice of the explicit form for the electric field 
is somewhat arbitrary, many other expressions having been proposed earlier (see e.g. Goldreich 
and Julian, 1969; Sutherland, 1979; Medin and Lai, 2008 and their references). We checked that 
the proton/ion saturation energy refiecting a balance between the acceleration and back reaction. 



and also the power of the different radiations are of the same order of magnitude when a simple 
electric quadrupole field is considered instead of the combination of the unsaturated and saturated 
fields as described above. 

The equations discussed above are valid only for protons. In case of accelerated ions, the 
relevant equation of motion can be arrived at by the replacements e — > Ze and m — > Am in 
(jlSp . We study in this section the case when the radiation intensity by the ions is the incoherent 
sum of individual nucleonic contributions provided that both photon and pion emission processes 
of the constituents become incoherent when the ions attain very large energies. Then one has to 
make the following modifications in the expression of V in (jlSp , 

• V-y — > ZV'y because only protons can emit photons; 

• Vt^o — > AVj^o because both proton and neutron emits vr'^ with the same probability and 
neutrons move together with the protons; 

• V.,^+ — > AVt^+ due to the radiation of vr^ by protons and the radiation of vr" by neutrons 
have same probability (from the point of view of neutrino radiation by charged pions we do 
not distinguish vr" and tt^). 

The equation of motion for heavy ions as detailed above is valid as long as the ion doesn't 
fragment. Fragmentation can occur as a result of the radiation recoil of the constituent nucleons. 
The numerical solution of the equation of motion shows that the transverse momentum transfer 
from the radiation k±_ ~ \^i^m exceeds the nucleon's binding energy (~ lOMeV for Fe ions) before 
the pion radiation occurs (see the right hand side of Figs. [4]). This means that the main part of the 
pion radiation arises from the acceleration of the individual proton constituents of the dissociated 
nucleus, as described by equation (fT5]) . Therefore, one doesn't need to deal with the details of ion 
fragmentation from the point of view of pion radiation, which would be given in this treatment by 
an incoherent sum of the radiation from the individual nucleon fragments. 

In the next subsection we shall present the numerical solution of the equation of motion for 
protons and the corresponding intensity of each type of radiations along the trajectory. 



3.2. Calculations and results 

One of the main goals of this calculation is to decide if one can find such values for the magnetic 
induction B and the rotation frequency P of the neutron star for which the power of charged pion 
radiation becomes comparable to the photon radiation. Therefore, based on the comparison of the 
timescales as presented in section 2, we focus our study on (S, P) pairs for which the intersections 
of the pion and photon cooling time curves with the acceleration time curve happen approximately 
at the same time. Based on the findings of the previous section (see Fig. [1]), we examined proton 
acceleration the following realistic magnetar configurations: (-B15, -P-3), (-B15, -P-2)) (-B14, P-s). 



The path of a proton is shown for (-B15, P-3) and different initial angles (0) measured relative 
to axis of the magnetic dipole on the left hand side of Fig. [2j It can be seen that the proton moves 
initially nearly along the magnetic field lines when the acceleration process begins, then shoves off 
the field lines at about z = 2-3 Rq and finally moves on a helical trajectory along the z-axis far 
away from the surface. 

On the top of the left hand side of Fig. [3l one can see that the maximum proton energy is 
almost independent of the initial angle when (-B15, -P-a)- This is also true in case of other values 
of B and P and also valid for the intensity maximum of the radiations as well. Therefore we 
only plotted the time dependence of the momenta and the power of radiations corresponding to 
d = 0.59pc (except the top left of Fig. [3|). 

It's interesting to observe time lags in the growth of both momentum components paral- 
lel/perpendicular (711 /7_l) to the field lines as can be seen on the top of Fig. [3] and on the right 
hand side of Fig. [21 in which the time dependence of the relative altitude ((r — Rq)/Ro) of pro- 
tons is shown. In this time region, the proton path is closer to circular (7_|_ 3> 711) and remains 
roughly constant because synchrotron radiation loss keeps pace with the transverse energy gain 
due to acceleration along the field. Then, the path changes to being more elongated due to the 
increasing electric field {\Eunsat\ ~ (r — Rq)/Rq), resulting in a change of the effective radius of 
curvature. This change of the path is indicated by a jump of the value of x = "7 /Re (at ~ 10"^'^ s) 
which can be seen in Fig. [U where we plotted the transverse momentum carried by the transfers of 
radiations {k± ~ mx^'^)- Between 10"^'' — 10~ s 7j_ is approximately constant and the effect from 
the proton energy gain is practically compensated by a decrease in the effective curvature radius, 
so the radiation rate is approximately constant (see the bottom parts of Fig. [3] and Fig. [3]). The 
increase in radiation rate after 10^^ s is due to the effective radius of curvature starting to decrease 
more slowly than the proton gains energy, until the saturated regime is reached where energy gains 
balance energy losses. 

From the bottom parts of Fig. [31 one can see that the vro radiation is irrelevant, moreover the 
charged pion radiation intensity becomes comparable to that of the photon radiation only in case 

of(i?i5,P-3). 

In Fig. [4l one can see that the transverse momentum transfer of the radiation is small com- 
paring to the total energy both for protons and for constituent protons of Fe ions. It supports 
the selfconsistency of our softness assumption for the radiation process implicit in our choice of 
the effective pion-nucleon interacton. However, k± is large enough to fragment the ions at the 
beginning of the trajectory as we noticed in the previous sections. 

The results of the numerical solution tell us that the protons can reach energies up to ~ 
10^/10^'^ GeV in the examined range of B and P. This is in good agreement with the results of the 
previous section where we compared the cooling time curves to the acceleration time curve. This 
conclusion is apparently not sensitive to the fact that the value of the momentum perpendicular 
to the magnetic field lines is almost the same order of magnitude as the parallel component due 
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to the essential presence of perpendicular components of the electric field. We also saw that for 
(Bi5, -P-3) or higher B and/or lower P it is the charged pion radiation which prompts the energy 
saturation of the protons (this is consistent with the conclusions obtained from the comparison of 
the characteristic timescales). 



4. Curvature Pion Neutrino Luminosity 

In terms of the Goldreich-Julian charged particle density ncj = (B/ecP) = 7x10^"^ Bi^P^^cm' 



^2 ^ 



for a pulsar or magnetar of radius a and period P s, with a polar cap solid angle fip = vri 
{27r'^a/cP) = 6.5 x lO^^agP""'^ sr, the outflow rate of ions with charge Z is 

Ni = a^%nGjZ-^c ~ 1.4 x lO^^^isa^Z^^P^^ g-i ^ (2I) 

For a vr"'" radiation efficiency r/jr at the ion terminal Lorentz factor, the resulting number of pion- 
decay t'^ and subsequent muon-decay v^ is 2 per proton, giving for a pulsar at distance D = 10 
Kpc a v^ or u^ flux at Earth of 

^y^c = ZNi/A-KD"^ ~ 1.2 X IQ-^'^'q^Bi^p-'^alD^^ cm^^g-i^ (22) 

where the label c denotes the curvature pion origin. The corresponding arrival rate in a detector 
of area A = IAq km^ is N^ = ^^A ~ 3.7 x lO^r/^Biso^p-^^o^ro^ km^^ yr^i. 

The typical energy of the muon neutrinos is determined by the energy losses incurred by the 
pions and muons before the decay. As the pions move out to radii comparable to the light cylinder, 
the vr^ Klein-Nishina energy losses are less important than the tt^ photon curvature losses, which 
occur on a timescale tc,7r ~ 1-2 x Id'^'^e^^a^P-^ s, where egvr = (evr/lO^GeV) is the pion energy 
normalized to 10^ GeV, while the vr^ decay time is tj^^T^ ~ 1.86x lO^eg^ s. The pions cool by curvature 
radiation until reaching the light cylinder at t^ ~ 1.7 x 10~^P_3S, where their energy has dropped 
to eg^ri ~ 1.9 X 10~ Og . Thereafter the pions move in the wind zone in essentially ballistic paths, 
cooling adiabatically on a timescale tad,-n = ^l{^txlI^-k)'^''^ = 1-4 x lO^^P^sCg eg^ s. The pions 
decay when their energy reaches eg^^ ~ 0.89 x 10~^P_3 Og ~ 10^ GeV. The corresponding Vp, has 
an energy e^,^ ~ e9,7r,d/3 ~ 33P_3 a^ TeV. The associated /u^ starts with (2/3) of the decaying pion 
energy and is subject to adiabatic losses. It decays on a timescale tp,d ~ 1-97 x 10^ eg^ s, where eg^ is 
the muon energy normalized to 10^ GeV, reaching at decay an energy eg^^ ~ 1.08 x 10~^P_3 Og ~ 
10^ GeV. The associated muon neutrino has an energy e^^ ~ (eg,^,d/3) ~ 3.3 TeV. 

The neutrino event rate in the detector, however, depends on the mixing of neutrino flavors 
because of oscillations while traveling from their source to Earth. Since the muon neutrinos from 
pion and muon decays are in two different energy regimes (33 and 3.3 TeV, respectively), one may 
consider their flavor oscillations independently. With a production flux ratio of fe:f^:^'T=0:l:0, for 
pion decay t'^, the observed ratio would be 0.2:0.4:0.4. For 3.3 TeV muon decay neutrino fluxes, 
with a production ratio of 1:0:0 and 0:1:0, respectively for v^ and P^, the respective observed 
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ratios would be 0.6:0.2:0.2 and 0.2:0.4:0.4. Since the neutrino detectors cannot distinguish between 
a muon neutrino and an anti-neutrino, the observed flux of 3.3 TeV muon neutrino would be a 
factor 0.6 times the production flux. The probability per v^ of resulting in an upward muon is 
Pu^fi ~ 1.3 X 10~^(ei/^/TeV), which for the 33 and 3.3 TeV neutrino cases lead to upward muon 
count rates from curvature pions, after taking into account oscillation effects, of 



^ _ f6.3 X lO^rt^Bi^pJ/^f^AoD^Q^e.JSSTeY) km-^yT-\ 
^'" ~ \9.5 X lQSnBi^PZl"'af'''A^Di;}{e,j2,.2,TeY) km-2yr-\ 



for an on-beam magnetar at 10 Kpc. 

In the case of proton acceleration, from Fig. [T]these would reach a maximum energy determined 
by 7r+ curvature losses for flelds i?^10^^ G and P ^ 10^^ s, or 515^^3 ~^- ^"^ these cases t^+<tc and 
the vr"'" radiation efficiency 7/,^ ~ 1, so the upward muon event rate should be detectable in the 1-100 
TeV range with cubic kilometer detectors, if the pulsars is on-beam and Bi^P_^ AqD^q >10~^, 
e.g. for sources in the local group. On the other hand, for values Bi^PZ^^ < 1, the pion efficiency is 
low, ry^ ~ tc/t-rr < 1, which would allow potentially detectable neutrino fluxes for sources in our own 
galaxy provided ^tt^IO"'^ — 10~^. For Bi^PZ^ <C 10~^, the steep behavior of the pion radiation 
efficiency implies an unobservably small number of events, even for galactic sources. Thus, there 
is a small range of magnetar field strengths i?i5>l and (short) periods P_3 ~ 1 for which these 
objects could be potentially detectable neutrino targets. 



5. Photomeson Neutrinos, Cosmic Rays and Photons 

Another mechanism for neutrino production in the polar caps is the photomeson process p7 — > 
vr"'" —> /^"''t'^ -^ ^'^i^e^fj.i acting on the magnetar surface X-rays, typically L^ ~ 10'^^L35 erg s~^ 
(Zhang et al, 2003). The photon number density near the surface is n^^^ = 1.8 x 10'^^L35^'^ 
cm~^, and for protons of energy much above the threshold e >> 0.1 — 0.2 GeV, the cross section 
(7p^ ~ 10"^*^ cni'^ implies a mean interaction time tp^ ~ 1.8 x 10^^L35~^'^ s, a mean free path 
tp^i ~ 5 X 10^L35~'^'^ cm, or an optical depth Tp^ ~ 0.18L35^/^a6. For L35>,1 and any but the 
fastest rotating pulsars, the photomeson interaction occurs on a length comparable or less than the 
light cylinder Ri = 5 x 10^P_3 cm, but much larger than the height h where acceleration balances 
pion or photon curvature losses. 

If the pion curvature process is inefficient, rj^^ <C 1, most of the protons survive beyond h, and 
inside the light cylinder they cool due to curvature radiation (or outside the light cylinder due to 
adiabatic cooling), until the photomeson process leads to a neutron and a vr"*" which takes 0.2 of 
the remaining proton energy. The pions and the subsequent muons then cool due to curvature (if 
inside Rl) and adiabatic (when outside Rl) losses in a similar way as in the previous section. For 
typical parameters, the average neutrino energy is < e^^ >~ 28P_3 L35~^'^^ag TeV, not much 
different from the value ~ 18P_3 Og TeV in the pion curvature initiated case of last section, since 
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both the protons (before undergoing pj) and the charged pions cool by photon curvature radiation 
inside the hght cyhnder, and adiabatically outside. The neutrino flux resulting from pj are larger, 
for Tp^ ~ 0.2 ^> t/tt, and since the typical neutrino energies from p'j are typically only ~ 2 larger 
than in the pion curvature case, in this 7/,^ <C 1 limit it would be hard to distinguish the pion 
curvature component from the dominant p'j flux. 

The situation is different for large !]■„ ^ 1. In this case it is mainly neutrons that propagate 
outwards beyond h. Being neutral, they do not suffer adiabatic losses, and undergo photopion 
interactions nj -^ it ~ on a timescale similar to protons, tny — 1.8 x 10^^L35^^'^ s. Typically this 
occurs inside or at the light cylinder. We take as a limiting example the case where this occurs 
near the light cylinder. The negative pions and muons then undergo adiabatic cooling outside 
Rl similarly to their positive counterparts, as discussed previously, resulting in P^. Ignoring the 
distinction between z/^ and f^, which cannot be discriminated in Cherenkov detectors, the resulting 
average neutrino energy is now < e^^ >~ 280L35~'^'^'^eg^^ TeV, where egni = (e„/10^GeV) is the 
initial neutron energy. The flux is ^^n-y ~ 5 x 10^^ PZ^ Bi^agLSb^'^ D^q cm~^s~^ around e,y ~ 280 
TeV from nj interactions, and the upward muon event rate, modulo the factors due to oscillation 
effect as previously discussed, is 

N^^n'y ^ 5.4 X 10^Si5Pr|a|L35^/'°Ao^ro'(e.^/280TeV) km-Vr-\ (24) 

for an on-beam magnetar at 10 Kpc, from wy interactions. The corresponding pion curvature 
neutrinos would lead to an upward muon event rate given by equation (f23l) . which is not too 
different. If observed, the difference in the average neutrino energies from pion curvature (~ 18 
TeV) and from wy interactions (~ 280 TeV) could help discriminate between the pion curvature 
and the photomeson production mechanisms. If there were such milisecond magnetars with 7/,^ ^ 1 
in our galaxy {D<W Kpc) this would imply, from eqs. (j23l [M|) upward muons event rates which 
would be strong in a completed cubic kilometer detector, and probably even in partially completed 
installations, so the current presence of such objects in our galaxy is questionable. 

These il, ■ B < caps will also result in a significant flux of escaping neutrons, which contribute 
to the cosmic ray flux. As an example, taking the case for the efficiency of vr^ radiation by protons 
from magnetars with -B15PZ3 = 1, for which r/j^ — t^+/tc ~ 0.2, the flux of neutrons at Earth and 
their energies are 

$„ ~ lO^^SisPrgDfo^ m"2 g-i SJ.-1 . e„ ~ 6 X lO^'^ eV . (25) 

A neutron of 6 x 10^'^ eV energy decays after traveling ~ 6 kpc. For comparison, the observed 
cosmic-ray number flux at this energy is ~ 9 x 10~^^ m^^ s^^ sr~^ (Nagano & Watson, 2000). 
However, the very small degree of observed cosmic ray anisotropy at these energies again suggests 
that the probability of finding such ry^r ^ 1 magnetars within our galaxy at any given time is very 
low. 

In the case of a newly born neutron star or magnetar, the escaping neutrons may interact with 
the expanding supernova-remnant (SNR) shell, as discussed by Razzaque, Meszaros & Waxman 
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(2003). For a hypernova which results in a niagnetar, the shell speed is f = W^^viq cm/s and 
reaches a radius Rsnr — 8.6 x lO^'^viot^ cm in t^ days. With a typical shell mass of Mgnr = 
2 X W^^msnr g, the column density of target atoms is Tja ~ 1.3 x lO'^^mgnrt^ v^q cm~^. The pp 
cross-section for ~ 10^^ eV energy incident neutrons is 100 mb and the corresponding optical depth 
is Tpp ~ lOmsnrt^ v^q . The energy of the neutrinos produced by a 7 x 10-'^'^ eV neutron via pp 
interactions ranges from ~ rriT^cP'^cm/'^ — 0.7 TeV up to ~ en/4 with a e,7^ distribution normalized 
to a multiplicity of 10~^ at this energy (Razzaque, Meszaros & Waxman, 2003). Here 7cm is the 
Lorentz factor of the center-of-mass of the pp interaction in the Lab frame. 



There may also be photo-hadronic interactions with photons in the SNR shell created from 
the SN explosion. The peak energy of these blackbody photons at creation is ~ 1^.2E^[ r-^^ keV 
for ~ 10^^£^5i erg SN explosion energy and ~ lO^^rn cm progenitor star's radius. In the SNR 
shell, however, they cool down to a peak energy of e^^sn — 2.I-E51 r^( v^q t^ eV. Thus neutrons of 
> 10^ eV energy may produce pions by interacting with them. Assuming that the total number 
of SN photons do not change after their creation, the column density of them in the SNR shell is 
S-y^sn — 3.65 X lO^^E'g^ r-^[ w{q t^ cm^ and the opacity for photo-hadronic interactions is < 1 with 
0.1 mb typical cross-sections at the resonances. The resulting neutrino energy in this case would be 
in the lO's of PeV range. There may be additional components of VHE ions accelerated in the SNR 
shell or the outer-gaps of a pulsar/magnetar which could also interact with the SNR shell (Guetta 
&; Granot, 2003). These interactions will produce additional neutrinos and hence identifying the 
neutrinos from pion curvature radiation may be difficult for t< few days after a magnetar is born. 

Ultra-high energy photons will also be produced, both through curvature radiation of charged 
pions and muons, and through decay of the associated neutral pions. E.g. pions produced by 
~ 10^ GeV protons would lead to pion-related UHE photons of luminosity Ljr^\ ~ lO^^Bi^PZ-^ Og 
erg s^^, decaying as the field drops and the period lengthens. The curvature photon energy is 
e^ ~ 3Pr3ag (evr/lO^GeV)^ PeV. These, as well as the muon curvature and neutral pion decay 
photons are all well above the 2meC^/ sin 9 threshold for one-photon pair production 'yB -^ e^ (e.g. 
Harding and Lai, 2006), where 6 is the angle of propagation relative to field direction. The optical 
depth above threshold for a path length lO^iJg cm is t^b ~ 5 x 10^°(S^f smOf/^Reie^/FeY)-^/^. 
(Another UHE photon opacity is photon splitting, a higher order mechanism which above threshold 
is less important than one-photon pair formation). The UHE photons will thus be degraded to 
energies below the threshold, e7<10 9zl MeV, for cap angles 9 ~ 10^^6'_i. 



6. Discussion 

Our discussion has assumed that the properties of the inner gaps near the polar caps of pulsars 
or magnetars allow protons to be accelerated to energies >6 x 10^ GeV. This is an open question, 
since the extent and properties of polar gaps (with Q ■ B < 0) including self-consistently pion 
radiation effects is unknown. Such studies, if undertaken, would have to take into account not 
only the effects of leptons acceleration from charged pion decays, but also high energy photons 
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from neutral pion decay and their subsequent cascades. Studies of space charge hmited gaps with 
electron acceleration and pair production (e.g. Harding and Muslimov, 2002; also Baring and 
Harding, 2002) suggest that Lorentz factors of order 10^ may be possible in fast-rotating, high field 
objects, in which case the effects discussed here may become important. 

The curvature pion radiation mechanism is likely to be of interest for fast rotating magnetars 
which accelerate protons. This is because, from Fig. [H we see that for magnetic field and period 
values i?i5P~3^1 the pion curvature radiation dominates the photon curvature radiative cooling. 
For such magnetic field and period combinations, the acceleration of heavy ions results in the 
fragmentation of the ions before pion emission occurs, which is expected only from the individual 
protons resulting from the fragmentation after they have been further accelerated. 

There is so far no observational evidence for millisecond periods among the handful of known 
magnetars (e.g. Kaspi, 2007; Woods and Thompson, 2004), most of which are in our galaxy and 
have periods of seconds. From eq. (j23p . we see that for i?>,10^^ G and millisecond periods, the 
muon event rate is so large that they might have been detected by now, if the magnetar were 
on-beam and at a distance D< Mpc. The likeliest explanation is that such millisecond magnetars 
do not exist in our galaxy (or else they might be off-beam). At the same time, it is reassuring that 
for periods P^l s the pion curvature efficiency predicted is rj-,^ = tc/^Tr '^ l, and the corresponding 
muon event rate in cubic kilometer detectors given by equation ()23p is negligible for known galactic 
magnetars. 

An interesting possibility is that some fraction of core collapse supernovae lead, at least ini- 
tially, to a millisecond magnetar. There is currently no observational evidence for this, but it 
is a plausible hypothesis, among others from fast convective overturn dynamo arguments, e.g. 
Thompson and Duncan (1996). In such objects, the envelope optical thickness to vN interac- 
tions is TyN ~ 10~^(ei,/10 TeV)(Me„t,/10MQ)(t;ent)/0.1c)~^(t/day)^^, and as long as the magnetic 
field and the rotation rate remain high, neutrinos produced by pion curvature and p7 interactions 
can escape without further reprocessing. Thus, core collapse supernovae which result initially in 
a millisecond magnetar with i?>10^^ G are expected to undergo a brief period of intense pion 
and ultra-high energy neutrino production, which significantly exceeds any electromagnetic energy 
losses. If the fraction of core collapse supernovae leading to such objects were 0.05 — 0.1, given the 
frequency of core collapse SNe in the LMC and in the local group, one could in principle expect 
some millisecond magnetars detectable in cubic kilometer detectors within timescales of years. In 
this scenario one would also expect (^ a significant <,10 MeV photon luminosity, which may be 
detectable by the GLAST GBM. 

In summary, we have pointed out that pions produced by protons interacting with the curved 
strong magnetic fields of magnetars may be an important energy loss mechanism for the accelerated 
particles, leading to secondary leptons and photons which can affect the properties of the inner 
gaps. Our model for the strong interaction processes involved here relies on the softness of the 
pion radiation process. Although the correctness of this assumption needs further, more detailed 



15 



investigation, this level of approximation is justifiable here, given the larger uncertainties in the 
astrophysical model. Our results indicate that the effects analysed in this paper could lead to 
copious neutrino production in the TeV-nPeV range, which would provide interesting targets for 
cubic kilometer scale detectors. 

This work is supported in part by NSF AST 0307376 and the Hungarian Science Foundation 
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at the Naval Research Laboratory. We are grateful to A.K. Harding and to the referee for useful 
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Fig. 1 . — Energy loss timescales t^+ and t^o for proton curvature emission of 7r+ and vr*^ , compared 
to photon curvature radiation foss timescale tc and acceleration time t^, as a function of proton 
energy (in GeV), for different field strengths B = lO^^Bis (G) and periods P = 10^'^P_3 (s). 
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Fig. 2. — On the left hand side: the proton's trajectory for initial angles relative the rotation axis 
O/Opc = 0.3, 0.4, 0.6 and i?i5,P_3 {9pc = 27.25°). 0pc is the angular size of the polar cap. On the 
right hand side: the time dependence of the relative altitude of protons above the surface of the 
neutron star for different values of B and P. 
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Fig. 3. — Acceleration of protons. The parallel and perpendicular momenta for i3i5,P„3 (a), and 
for Pi5,P_2 (^pc = 8.33°) and Bu,P_3 {9pc = 27.25°) where 6 = 0.5^ (b). The intensity of each 
type of radiation for i?i5,P_3 (c), for i?i4,P_3 and S15, P_2 (d) where 6 
lines indicate the switch between the unsaturated and saturated fields. 
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Fig. 4. — Transverse momentum transfer the radiations from protons (left hand side) and con- 
stituent protons of Fe ions (right hand side) along the path. The photon radiation dominates in 
case of (i?i5, -P-2) and (-B14, i^-3), and in case of (-B15, -P-3) up to t ~ lO^^s. 



